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ABSTRACT 



Context. Recent spectro-polarimetric observations have provided detailed measurements of magnetic field, velocity and 
intensity during events of magnetic field intensification in the solar photosphere. 

Aims. By comparing with synthetic observations derived from MHD simulations, we aim to discern the physical pro- 
cesses underlying the observations, as well as to verify the simulations and the interpretation of the observations. 
Methods. We consider the temporal evolution of the relevant physical quantities for three cases of magnetic field inten- 
sification in a numerical simulation. In order to compare with observations, we calculate Stokes profiles and take into 
account the spectral and spatial resolution of the spectropolarimeter (SP) on board Hinode. We determine the evolution 
of the intensity, magnetic flux density and zero-crossing velocity derived from the synthetic Stokes parameters, using 
the same methods as applied to the Hinode/SP observations to derive magnetic field and velocity information from the 
spectro-polarimetric data. 

Results. The three events considered show a similar evolution: advection of magnetic flux to a granular vertex, devel- 
opment of a strong downflow, evacuation of the magnetic feature, increase of the field strength and the appearance of 
the bright point. The magnetic features formed have diameters of 0.1-0.2^^ The downflow velocities reach maximum 
values of 5-10 km/s at r = 1. In the largest feature, the downflow reaches supersonic speed in the lower photosphere. In 
the same case, a supersonic upflow develops approximately 200 s after the formation of the flux concentration. We find 
that synthetic and real observations are qualitatively consistent and, for one of the cases considered, agree very well 
also quantitatively. The effect of finite resolution (spatial smearing) is most pronounced in the case of small features, 
for which the synthetic Hinode/SP observations miss the bright point formation and also the high- velocity downflows 
during the formation of the smaller magnetic features. 

Conclusions. The observed events are consistent with the process of field intensification by flux advection, radiative 
cooling, and evacuation by strong downflow found in MHD simulations. The quantitative agreement of synthetic and 
real observations indicates the validity of both the simulations and the interpretations of the spectro-polarimetric 
observations. 
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1. Introduction sure deficit, the gas inside the flux tube will then be pressed 

^ together (together with the frozen-in magnetic field) by 

Magnet ic field is ubiquitou sly present in the solar photo- the surrounding gas, causing the magnetic pressure to in- 

' sphere (|de Wiin et alJl2008D . On granular scales, it under- crease until a balance of total pressure (magnetic + gas) is 

goes continual deformation and displacement. It is swept by reached. 

the horizontal flows and concentrated in the intergranular , . „ - ^ - u ^ ^ 
, , , , . 1 n 1 1 . 1 . . 1 -L fie convective collapse extents the concept to the con- 
lanes. I^ lows are able to compress the neld so that the mag- . , T-i-, . . -.1 n ^ u • ^1 1 
, 7-,9 /o 1 .1 1 • vective mstability. It starts with a nux tube m thermal 
netic energy density B^/S tt approaches the kmetic energy -, u • i -i-u • -^u u j 
J 2 /o £ ^if a H t^ ^ J l mr-o l hur • \ Unrr [\ rru- mcchauical equilibrium with the surrounding hydro- 
density pv^/2 of the flow (IParkerlll963l : IWeisslll966l ). This , ^ . i . --n -> i c- / i 

. ' n ~: — — statically superadiabaticly stratmed plasma. Since external 

results m a magnetic neld strength ol a lew hundred Gauss . ^-n \- • -- i ^ui --i ^• 

, , p .1 . , .n . 1^ . .1 stratmcation IS convectively unstable, any vertical motion 

at the solar surface, l^urther intensification to kG strength -^u • a ^ u u t-c j T-k ^ a -n 

1 . 1 ^1 1 . r 1 X T 1^- withm the flux tube can be amplified. Downward flow will 



is driv en by t he mechanism referred to as: s uyeradiahatic . i-, -, jj-fL ^-ir a 

^ \ ^ / — irr^rroN ^- 7 7 i o t-. i ^ grow m amplitude and dram the material Irom the nux 

efject (IFarkerl 1978i ), convective collapse (I Webb Roberts. ? ^ m^ ^- • -xi 



i V^r^o l In r7 1 — qi -, J — TT'^ tube. Ihe process continues until a new equilibrium with 

119781 : Ibpruit fc Zwei bell 119791 ) or convective intensification , n • i j T\-ai . ^ r . i ^ 

' ^ ^1 X ilhJAoK a strong field is reached. Diflerent aspects ol th e concept 

(IGrossmann-Doerth et al.lll998l ). i u-^r^ • Ic u i J 

^ ^ have been the subj ect ol extensive research (see Ibchiisslerl 

The first two concepts are a theoretical idealization of 1 19901 ISteinerl [l999l for reviews). The term convective in- 
the process. The s uperadiabatic effect contains the basic tensification is used for magnetic field intensification in re- 
idea. iParkerl (|l978l ) pointed out that a thermally isolated alistic MHD simulations, where the process occurs in its 
dowfiowing gas within the flux tube in a superadiabatically full complexity. It is driven by the thermal effect in the 
stratified environment will be accelerated, which would lead surface layer of the magnetic concentration. There, due to 
to evacuation of the flux tube. Because of the resulting pres- the presence of magnetic field, heat transport by convection 
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Fig. 1. Maps of the whole simulation domain at t = 140 s. Normalized continuum intensity at 630 nm(left), vertical 
component of magnetic field (middle) and velocity (right) at a geometrical height roughly corresponding to the level of 
('?"50o) = 1 (^ 930 km above the bottom of the computational box). 



is reduced. The material inside the concentration radiates 
more that it receives. This leads to cooling of material which 
thus starts to sink and partial evacuation of the concentra- 
tion occurs. Contraction of the magnetic concentration by 
the surroundings (result of the pressure imbalance) leads 
to an inc rease in n i agnet i c field strength. Thus, the sim- 
ulations (iNordlundl Il983l: IG rossmann-Doerth et al.l 119981 : 



iGadun et all 120011 : IVogler et al.ll2005l : ICheung et ^l200? 

bear out the basic properties described by idealized con- 
cepts. That is the downflow, the evacuation of the magnetic 
structure, the field increase and, in some cases, establish- 
ment of a new equilibrium. The 3D MHD simulations show 
that the strong magnetic concentrations form as the hori- 
zontal flows in the intergranular lanes advect weak, nearly 
vertical field and concentrate it at the vertices of granu- 
lar and m esogranular downflow lanes (|Stein fc Nordlundl 
I1998LI2006D . Larger magnetic structures form at sites where 
a granule submerges and the surrounding field is pushed 
into the resulting dark region. Whether the formed concen- 
tration appears dark or bright in the continuum intensity 
depends on whether the vertical cooling is compensated 
or not by the la teral heating d ue to horizont al energy ex- 
change (|Bercik et a l. 2003: Vo gler et al.l[2005l ). This forma- 
tio n scena r io is c onsis tent with the obse r vation s described 
bv iMulled (|l983r ) and iMuller fc Roudied (|1992D . Their ob- 
servations show that network bright points form in inter- 
granular spaces, at the junction of converging granules as 
the magnetic field gets compressed by the converging gran- 
ular flow. 



The 2D simulations by lGrossmann-Doerth et"al] (|l998l ) 
revealed that magnetic flux concentrations formed by con- 
vective intensification can evolve in different ways. They 
present two possible outcomes. Depending on the ini- 
tial magnetic flux, a magnetic concentration can reach a 
stable state after the process, or can be dispersed due 
to an upflow that develops as high speed downflowing 
material rebounces from the dense bottom o f the tube. 
Similar results were prese nted by iTakeuchl (|l999l ) and 
ISheminova fc GadunI (|2000f ). 



Observational ev i dence was found for both cases. 
iBello Gonzalez et al.l (|2008l ) reported on the formation of 
a magnetic feature at the junction of i ntergranular lanes, 
without any significant upflow observed. iBellot Rubio et al.l 
(2001) , on the other hand, detected a strongly blueshifted 
Stokes V profile originating in a upward propagating 
shock, 13 minutes after the ani p lificat ion of magnetic field. 
ISocas-Navarro fc Manso Saind (|2005l ) found that super- 
sonic upflows are actually quite common. 

Events that are interpreted as convective c ollapse were 
detec ted also with the spectropolarim eter fSP) (|Lites et al] 
l200lf ) of the Sola r Optical telescope (iTsuneta et al.ll2008D 
on bo ard Hin ode ()Kosugi et al. 20071 ). Both. Shimizu et al.l 
(|2008D and iNagata et all (|2008D show cases of high 
speed downfiows followed by magnetic field intensifica- 
tion and bright poin t appearance. The event described by 
^Na gata et al.l (|2008l ) shows stronger field strength and up- 
flow in the final phase of evolution. 

In this paper we give three examples of magnetic field 
intensification from MURaM simulati ons and make de tailed 
comparison w i th th e results of .Nagata et al.l (|2008f ) and 
IShimizu et all (|2008h . 



2. Simulation data and spectral synthesis 

We use 3D radiative MHD simulations of a thin layer con- 
taining the so l ar surface carried ou t with the MURaM code 
(|VQglerl [2OO3I : IVogler et all l2005f ) in a 6 x 6 x 1.68 Mm 
domain with non-grey radiative transfer included. Vertical 
and horizontal spatial resolution is 10 km and 14 km re- 
spectively. The bottom and top boundaries are open, per- 
mitting free in and outflow of matter. The initial magnetic 
field of {\B\) = 200 G, is introduced in a checkerboard-like 
2x2 pattern, with opposite polarity in adjacent parts. As 
the field is redistributed by convective motions, opposite 
polarities are pushed together and dissipated, so that the 
unsigned magnetic flux decreases with time. In this way, 
the run simulates the decay of the mag netic field in a mixed 
polarity region. Local dynamo action ([Vogler fc Schiisslerl 
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Fig. 2. Evolution of the 
continuum intensity at 
630 nm and the magnetic 
field in a 4" x 3'' sized 
region. Left-hand side 
(double column): original 
spatial resolution; red and 
blue contours outline the 
downflow of 6 km/s at 
80 km and upflow of 4 km/s 
at 400 km above r = 1, 
respectively; horizontal 
lines mark the positions 
of vertical cuts shown in 
Figs. [4] and [H boxes corre- 
sponding to cases I, II and 
III at coordinates [2" , 1.5''], 
[0.5'',!''] and [3.5'',!''], 
respectively, enlarged in 
Fig. [3l Right-hand side: 
synthetic Hinode/SP obser- 
vations; left column shows 
continuum intensity while 
the right column shows 
the apparent field strength 
(see the text); red contour 
marks the region with 
0.01 pm of signal excess 
(see text); yellow crosses 
mark the positions of pixels 
that we studied in detail in 
Fig.O Grey scales cover the 
range of 0.6-1.6 and 0.8-1.2 
for normalized intensity 
at original and reduced 
resolution, respectively, 
±1000 G for the vertical 
component of magnetic 
field (original resolution) 
and ±200 Mx/cm^ for the 
apparent longitudinal mag- 
netic flux density (Hinode 
resolution). 



I2007D does not occur since the magnetic Reynolds number 
is below the threshold for dynamo action. We do not expect 
that the intensification process would be significantly dif- 
ferent at higher Reynolds numbers (smaller grid cells). The 
simulated magnetic flux concentrations considered here are 
well resolved. The main effect of an increased resolution 
would be a decrease in the width of the boundary layers 
between the flux concentration and the surrounding down- 
fiows. 

We examine a 30 min sequence of simulation snapshots 
with a cadence of approximately 90 s. The snapshots shown 
in this paper have {\B\) ^ 150 G at r = 1. In Fig. (TJ we 
show the continuum image and maps of the vertical com- 
ponents of the magnetic field and velocity for a snapshot 
at time t = 140 s (t = s corresponds to the first snapshot 
considered in this paper). 

In order to synthesize the Stokes profiles, the physical 
parameters from the simulation are used as a n input for the 
ID LTE radiative transfer code, SPINOR (jFrutiger et all 



I2OOOD . A spectral range that contains Fe I lines 630.15 and 
630.25 nm is sampled in steps of 7.5 mA. The Fe abun- 
dance used for the synthesis has been taken from iTheveninI 
(1989) and the v alues of the oscillator strengths from the 
VALD database (iPiskunov et al.lll995f ). Before comparing 
with Hinode/SP (|Lites et al.l 200l[ ) observations, the syn- 
thetic line profiles have been treated to bring them to the 
same resolution as Hinode data. Firstly, a realistic point 
spread function PSF (^Danilovic et al.ll2008f ) has been ap- 
plied to the synthesized Stokes profiles. The PSF takes 
into account the basic optical properties of the Hinode 
SOT/SP system and a slight defocus which brings the con- 
tinuum contrast of the simulation to the observed value of 
7.5%. Secondly, to take into account spectral resolution of 
the spectropolarimeter, the profiles are convolved with a 
Gaussian function of 25 mA FWHM and resampled to a 
wavelength spacing of 21.5 mA. Thirdly, noise of 10 Ir 
is added. Velocities are determined from the shift between 
the Stokes V zero-crossing wavelength and the line core 
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Fig. 3. Enlargement of the regions, having a size of 0.75^' x 
0.75^', outlined by yellow squares in Fig. O From left to 
right: cases I, II and III. Continuum images (left column 
for each case) and maps of the vertical components of the 
magnetic field (right) with overplotted horizontal velocities 
at approximately 80 km above r = 1 (white arrows). Red 
contours outline locations of dowflows of 5 km/s or larger 
and yellow contours the vertical component of vorticity ex- 
ceeding ±0.2s~^ (cases I and II) or ±0.1s~^ (case III). 

position of the Stokes / profile averaged over pixels with 
polarization signal amplitudes less than 10~^/c- We cor- 
rect these velocities for convective blu eshift by subtract ings; 
150 m s~^. Finally, the procedure bv lLites et al.l (|2008l ) is 
used to calculate the longitudinal apparent magnetic flux 
density (shown in Fig. [2]). 

3. Results 

Fig. [2] shows the evolution of the continuum intensity and 
the magnetic field in a 2.9 x 2.2 Mm subdomain of the 
simulation, over approximately 7 minutes. Shown are maps 
at the original resolution of the simulation (left-hand side) 
and at the spatial resolution of Hinode (right-hand side). 
For the original resolution, we show the vertical compo- 
nent of the magnetic field from the simulations (right col- 
umn on the left-hand side ). The longitudinal apparent flux 
density (|Lites et al.ll2008l ) retrieved from synthetic Stokes 
profiles is shown at the Hinode resolution (right column on 
the right-hand side). During the period shown, three bright 
points appear near the coordinates [2'\ 1.5^'], [0.5^', 1^'] and 
[3.5'', 1''], in the regions outlined by yellow squares. We re- 
fer to them here as cases I, II and III, respectively. They 
are identified in the top left frame. 

3.1. Horizontal flows 

In all three cases, magnetic field is advected by the flow to 
the junction of multiple granules, where it is confined and 
concentrated. The evolution of magnetic concentrations is 



examined more closely in Fig. [3l where the regions out- 
lined by yellow squares in Fig. O are enlarged. Horizontal 
velocities, at approximately 80 km above (r) = 1, are rep- 
resented by arrows in continuum maps and magnetograms. 
Regions with a strong vertical component of vorticity are 
also outlined. 

Vortex flows around strong downflows are quite com- 
mon in simulations of granulation (yNordlund il986 ) . They 
are formed at the vertices between multiple granules, where 
flows converge and angular velocities with respect to the 
center of downflow increase due to angular momentum con- 
servation. The lifetime of vortex flows depends on the dy- 
namical behavior of the neighboring granules. In a 30 min 
run that we examined, they last from less then 90 s to 
approximately 10 min. The ones shown in Fig. [3] are short- 
lived. In case II, the swirling motion persists for at least 150 
s, until the shape of the neighboring granules is changed. 
The formed flux concentration is then squeezed between 
two granules and stretched into a flux- sheet- like feature. 
Later on, as the granules evolve, the magnetic feature is ad- 
vected to the left and trapped again in the junction between 
newly formed granules. In case HI, a strong preexisting flux 
sheet is carried by the flow toward the junction of the gran- 
ules. There, the field is caught in a vortex flow that lasts 
for at least 100 s. The flow becomes disturbed and then 
starts to swirl in the opposite direction. The vortex axes 
are inclined, in both cases II and HI, with respect to the 
vertical direction. Also in both cases, there is a spatial and 
temporal coincidence between the existence of vortex flows 
and the formation of high speed downflows near the surface 
layer. 

The vortices shown in Fig. [3] are less then 0.25'' in di- 
ameter, much sma ller then the ones found in observations 
(B onet et al]l2008[ ). The size of the observed vortex flows is 
of order l"and the average lifetime is around 5 min. They 
seem to outline supergranulation and mesogranulation cells. 
As the authors of the observational study suggest, most vor- 
tices could have been missed because of the limited spatial 
resolution (Swedish Solar Telescope) and the method used 
(tracking the motion of bright points in G band). A short- 
lived vortex flow preceding formation of a netw ork bright 
point has been observed by iRoudier et al.l (|l997l ) . 

No vortex is detected during the formation of flux con- 
centration in the case I. In this case, a strong downflow, and 
later upflow develops inside the flux concentration, whose 
shape and contrast in the intensity map is greatly affected 
by the convective motions. Spatial and temporal fluctua- 
tions across the magnetic feature make it difficult to define 
its center. The bright points, visible in continuum maps, 
form where fragments of flux concentrations become nar- 
row so that lateral heating becomes important and in the 
regions with higher field strength, i.e. higher evacuation. 

During their further evolution after the intensification 
phase {t > 420 s), the flux concentrations in cases I, II and 
HI have different fates. In case II, it meets a magnetic fea- 
ture of opposite polarity and vanishes in ~ 3 minutes. The 
feature formed in case I is gradually stretched by the flow 
and fragmented after approximately 9 minutes. A few min- 
utes after that, the larger fragment encounters the opposite- 
polarity flux concentration from case HI. The cancellation 
continues for ~ 12 minutes until both features disappear. 
The smaller fragment of the feature from case I is caught in 
a granular vortex flow approximately 5 minutes later which 
leads to further intensification. 
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Fig. 4. Time sequence of 
vertical cuts through the 
magnetic concentration cor- 
responding to case I (posi- 
tions of the cuts are marked 
by yellow horizontal lines 
in right-hand images in 
Fig. [2). The plotted cuts 
correspond to the same in- 
stances of time as sam- 
pled in Fig. O From left 
to right: density excess with 
respect to the the mean 
density at every height, ver- 
tical component of mag- 
netic field, vertical compo- 
nent of velocity (overplot- 
ted arrows show the com- 
ponents of the velocity field 
in the x-z plane), and radia- 
tive heating rate. Horizontal 
solid lines follow the levels 
of log r = 0,-1 and —2. 
Vertical dashed lines mark 
the position of the height 
profiles given in Fig. [5l 



3.2. Field intensification 

The evolution of the magnetic structures formed during the 
studied time interval can be followed in more detail by tak- 
ing vertical slices through the regions marked by horizontal 
lines in Fig. [3l We will examine three cases. 

3.2.1. Case I 

Vertical cuts through the regions corresponding to this case 
are given in Fig.m Shown are the density excess (relative to 
the mean density at each geometrical height), the vertical 
components of magnetic field and velocity, and the radia- 
tive heating rate. Cuts at the instant t = show the fluc- 
tuation in density outlining the granulation structure, with 
less dense, hotter granules and a cooler intergranular lane, 
where magnetic field of a few hundred G is accumulated and 
a downflow of a few km/s is present. The fluctuation in ra- 
diative heating rate is highest in the thin layer near optical 
depth unity, where upflowing granular material cools most 
strongly. The regions along the intergranular lane seems 
to undergo cooling in the layers above r = 1. In the next 
instant, 140 s later, a strong magnetic concentration has 
already formed. Variations in magnetic field strength at 
a given geometrical height are accompanied by variations 
in the gas pressure. Regions with increased field strength 
have lower density, which, in turn, shifts the level of opti- 



cal depth unity downwards relative to neighboring regions 
with weaker field. The depression of the visible surface in- 
side the magnetic structure gives rise to radiative heating 
through the sidewalls from the hot neighboring material. 
The downflow velocity is significantly increased from t = 
to t = 140 s, becoming supersonic in the layer between 100 
km above and 350 km below the surface. By t = 280 s it 
is greatly diminished again inside the flux concentration. 
The diameter and the shape of the flux concentration con- 
tinuously evolves and does not reach a stationary state. 
At t = 280 s strong radiative heating is present inside the 
concentration. This is also the moment when continuum 
intensity is the highest as shown later (Fig. [9|). 

In the last two sets of snapshots, shown in Fig. [H an 
upflow is visible inside the flux concentration, extending al- 
most from the bottom of the simulation domain. It reaches 
supersonic velocities in the upper layers of the photosphere. 
Its direction is inclined with respect to the vertical. The 
upflowing material inside the concentration exhibits strong 
cooling at the surface. As the material refills the region in- 
side the flux concentration, the level of optical depth unity 
shifts upwards and the magnetic field strength decreases. 
However, the upflow does not lead to a complete dispersal 
of the field, as described in the last paragraph of Sec. 13. li 

The evolution of this flux concentration is further illus- 
trated by the vertical profiles given in Fig. [5l Three phases 
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Fig. 5. Height profiles of various quantities in the flux con- 
centration corresponding to the case I (magnetic field, ver- 
tical component of velocity, density and temperature) at 
the positions marked by dashed lines in Fig. [4] at different 
time instances. Thick vertical lines mark the positions of 
the optical depth rsoo = 1 for each atmosphere. Thin ver- 
tical lines in the top two panels mark the positions of the 
height of formation of Fe I 630.25 nm line (see text). 



before and after the field intensification are shown. The 
horizontal positions of the profiles are marked by verti- 
cal dashed lines in Fig. [H The profiles of magnetic field 
strength, vertical component of velocity, density and tem- 
perature for each phase are shown, together with the cor- 
responding level of optical depth unity, marked by thick 
vertical lines. The level h = corresponds to the level of 
the mean optical depth unity (rsoo) = 1 for the whole snap- 
shot. Thin vertical lines in the upper plots mark the posi- 
tion of the average height of formation of the Fe I 630.25 
nm line, which is defined by calculating the centroid of 
the contribu tion function (CF) for line intensity depression 
([Solanki fc Bruls 1994), at the wavelength of the line core. 
The definition gives only a rough estimate of the heights 
sampled, since the lines are formed over a large portion of 
the photosphere and the CFs are asymmetric in wavelength 
due to the presence of high velocity gradients. 

The figure shows that, in the interval from t = s to 
t = 140 s, the downflow extends to deeper layers and in- 
creases in amplitude, reaching 10 km/s at 200 km below 
the surface. The magnetic field strength increases from a 
few hundred to more than 2000 G at the level of rsoo = 1- 
The significant reduction in density results in a shift of 200 
km in the optical depth unity level. As a result of radiative 
heating the temperature gradient is flatter and the temper- 
ature is around 400 K higher at = 0. The temperature 
difference at equal optical depth is even larger, exceeding 
1000 K at Tsoo = 1- The estimated height of formation of 
Fe I 6302 shifts from roughly 300 km to 150 km above the 
surface when mag netic field concentration forms , which is 
in agreement with iKhomenko fc ColladosI (|200 7l) . This ef- 
fect contributes to the increase of observed magnetic field 
strength in the first two phases. 




h [km] h [km] 

Fig. 7. The same as Fig. O but for case II. 




-200 200 400 600 -200 200 400 600 

h [km] h [km] 

Fig. 8. The same as Fig. [5l but for case III. 



In the third phase, the upflow leads to a density en- 
hancement and a weakening of the magnetic field by a 
few hundred Gauss in the lower photosphere. The level of 
Tsoo = 1 shifts upward by 100 km. No discontinuity is visi- 
ble in the vertica l component of velocity, in contr ast to the 
event studied bv lGrossmann-Doerth et all (|l998l ). 

3.2.2. Cases II and III 

Vertical cuts through the regions corresponding to these 
cases are shown in Fig. [6l The positions of the cuts are 
marked by horizontal lines in Fig.[3l The parameters shown 
are the same as in Fig. IH 

At t = 0, weak magnetic structures in the intergranular 
lanes are already present in both cases. Strong downflows 
along the magnetic field lines persist over the entire time 
interval shown. They reach supersonic velocities in some in- 
stances, but only in the layers below r = 1. In case II, the 
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Fig. 6. The same as Fig. |4] 
but for cases II and III. 
Each frame is subdivided 
into a left part displaying 
case II and a right part 
showing case III. Vertical 
lines mark the positions of 
the height profiles given in 
Fig. [71 and Fig. El 



downfiow is strongest at t = 140 s, just before the vortex 
flow becomes deformed owing to evolution of the neighbor- 
ing granule. Then it gets somewhat reduced at t=280s, as 
the structure is squeezed between the granules. Later, at 
t = 380 s, a strong downfiow is present again in the newly 
formed granular junction, to which the flux concentration 
is carried and where it is confined again. The rightmost 
columns show strong cooling of the layers above the level 
of '?"500 = 1 in both cases. At t = 420 s, the interior of the 
flux concentrations reaches radiative equilibrium. At that 
moment, the flux concentrations in both cases are almost 
vertical, with diameters of approximately 0.1^' at rsoo = 1- 

Figures [71 and [8] show height profiles along the positions 
marked by the dashed lines in Fig. [6l Since the flux concen- 
trations are inclined, the choice of curved lines, along the 
axes of magnetic structures, would be more logical here. 
However, we are interested in what would be observed at 
the disc center and the observables (i.e. the Stokes profiles) 
depend on the vertical profiles of the physical parameters 
which are shown here. In the case II, we choose to show 
profiles at instants t = 0, 280 and t = 420 s (solid, dotted, 
dashed respectively). For the case III, profiles at the last 
three instants from Fig. [6l are shown. Both cases exhibit 
a similar behavior: there is a significant evacuation of the 
magnetic concentration due to the strong downfiow and an 
increase of the magnetic field strength from a few hundred 



G to kG values in the process. There is a persistent dowfiow 
that reaches 5 km/s (in the vertical direction) at rsoo = 1 
and a decrease in density with a corresponding shift in the 
height of Tsoo = 1. In both cases, the shift in rsoo = 1 
reaches 150 km and the magnetic field strength increases 
to 2000 G at that depth. The height of formation of the 
Fe I 630.25 nm line shifts to deeper layers as the magnetic 
field strength increases. In the last instant, the Fe I 630.25 
nm line probes layers with kG fields and vertical velocity 
of 5 km/s, in both cases. 



3.3. Comparison with synthetic Hi node observations 

The intensity contrast is considerably reduced at Hinode 
spatial resolution, as can be seen in Fig. [2l (note that the 
color scales are adjusted). The maps of longitudinal appar- 
ent magnetic flux density are shown for comparison with 
the magnetograms at the original resolution. Fine struc- 
ture in both maps, continuum and magnetograms, is lost 
due to smearing. The bright points that correspond to cases 
II and III are barely discernible. 



3.3.1. Case I 



In ord er to compare case I with the results of lShimizu et al.l 
(|2008[ ). we calculate the parameter they refer to as signal 
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Fig. 9. Temporal change of normalized continuum intensity (left), vertical component of magnetic field (middle) and 
vertical component of velocity (right), for cases I (top row), II (middle row) and III (bottom row) at original (stars/solid 
line) and Hinode resolution (crosses/dashed line). Vertical component of magnetic field and velocity at log rsoo = — 2 (for 
original resolution) are each plotted together with the longitudinal apparent magnetic flux density (note the separate scale 
employed for in the middle column) and velocity retrieved from Stokes V profiles (Hinode resolution), respectively. 






Fig. 10. Synthetic Stokes I (left) and V (right) profiles at 
original (upper row) and Hinode resolution (lower row) for 
the times t = 0, 140 and 420 s marked by dotted, solid 
and dashed, respectively. The intensity is normalized to the 
mean continuum intensity in the full snapshot domain. 



excess (SE). It is defined as the Stokes V profile integrated 
over the spectral range of 250 — 400 mA redward from the 



nominal line center. This wavelength range corresponds to 
dowfiows of 7 — 14 km/s, so that high values of SE indicate 
the presence of hight-speed downflows. The red contour in 
the right-hand side of Fig. [2] outlines a region with SE of 
0.01 pm, which is a factor of 10 higher than in the surround- 
ings. It fits the location and time of the strong downflows 
visible at the original resolution (left hand-side of Fig. [2l). 
Both, synthesized (here) and observed (Shimizu et al. 2008) 
SE occur simultaneously with the appearance of the bright 
point and the intensific ation of the magnetic field. However, 
the event described bv IShimizu et al.l (|2008l ) (Fig. 8) lasts 
at least 6 minutes, while in our simulations it is present for 
less than 4 minutes. 

Figure [9] compares the observable parameters with the 
corresponding values at original resolution. The locations 
chosen for the plots are marked by yellow crosses in Fig. [2l 
The pixels are selected such that regions of downflow and 
upflow are covered, as well as the evolution of the bright 
point. The plots show the temporal change of normalized 
intensity, longitudinal apparent flux density (|Lites et al.l 
2008), and zero-crossing velocity retrieved from Stokes V 
profiles of the Fe I 630. 25 nm line, so t hat a direct compar- 
ison with the results of iNagata et al.l (l2008h can be made. 



Overplotted (stars/solid lines) are the normalized intensity 
at the original resolution, the vertical component of mag- 
netic field and velocity at logrsoo = — 2. 

The comparison shows that the correlation between the 
parameters at original and reduced resolution is not so 
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good. This is a consequence of different effects. Firstly, 
since the formation region of the Fe I 630 nm lines ex- 
tends over a large part of the atmosphere, there is smearing 
along the line of sight and the resulting profiles carry in- 
formation from different layers with different values of the 
physical parameters. Secondly, the smearing due to spatial 
resolution affects every pixel differently depending on its 
neighboring pixels. The profiles could be more asymmetric 
and broader due to averaging over inhomogeneous regions 
and the amplitude of Stokes V can be affected. Figure [TOl 
shows examples of Stokes / and V profile at original and 
reduced resolution at t = 0, 140 and 420 s. The profiles 
show stronger asymmetries and higher Doppler shifts at 
original resolution. The amplitudes of Stokes V profiles at 
Hinode resolution are reduced by a factor of 10. They are 
of t he same order of ma gnitude as the observed ones (Fig. 
2 of iNagata et al.l (|2008l )). Both, synthesized and observed 
Stokes / profiles show an extended red wing at the phase 
with strong downflow. 

Figure [9] reveals that intensity contrast is significantly 
reduced by spatial smearing. Thus at t = 230 s the contrast 
of the bright point drops from 100% at original resolution 
to below 30% at the reduced resolution. The velocity de- 
termined from Stokes V profiles is at most times somewhat 
lower than the velocities at logrsoo = —2. However, the 
general trends of the paran ieters are preserved . For both, 
simulation and observation (|Nagata et al.ll2008l ). the inten- 
sity reaches its maximum approximately 100 s later than 
velocity. This is consistent with the idea that the enhanced 
brightness results from the evacuation due to the down- 
flow. When the velocity changes sign, the brightness starts 
to fade. The maximum brightness of the simulated and ob- 
served bright points are similar after degrading the sim- 
ulated intensity to Hinode resolution, while the velocities 
derived from the simulations are consistent, but somewhat 
smaller. They reach approximately 5 km/s (downflow) and 
0.5 km/s (upflow), while values up to 6 km/s and 2 km/s, 
respectiv ely, are observed. Th e magnetic field strength ob- 
tained bv lNagata et al.l (|2008l ), by inversions, increases from 
a few hundred to 2000 G in approximately 200 s, which is 
in agreement with simulations. 

3.3.2. Cases II and III 

The middle and bottom rows in Figure [9] show the tempo- 
ral evolution for cases H and HI, respectively. We find a 
simultaneous increase in all parameters (continuum inten- 
sity, apparent flux density and zero-crossing velocity). In 
both cases, the downflow is suppressed at t = 280 s pos- 
sible because of the evolution of the surrounding granules. 
The velocity at logrsoo = —2 reaches 6 and 4 km/s for 
case II and HI, respectively. As the magnetic field strength 
increases, the intensity follows and, in both cases, the in- 
tensity contrast reaches 50% with respect to the mean at 
original resolution. 

The influence of spatial smearing is on the whole similar, 
although a few differences relative to case I are noteworthy. 
Firstly, due to the small size of the features, the correlation 
of the magnetic and brightness signals is largely destroyed. 
Although the magnetic signals are comparable to the values 
in case I, the brightness is lower, so that the features are in- 
conspicuous in continuum images. The increase in magnetic 
flux density is followed by just a small fluctuation in inten- 
sity. Secondly, although the magnetic field strengths and 



diameters of the features are comparable for cases II and 
HI at original resolution, the case II shows lower magnetic 
flux density after convolution. Thirdly, the signature of the 
strong dowflow in case II is lost due to spatial smearing. 



4. Summary 

Our case study of three examples of magnetic field intensifi- 
cation has shown that in all three cases, the field is advected 
to the junction of several granules. There, it is confined 
by converging granular flows, which, in two cases, form a 
vortex. Owing to the presence of the magnetic field, the 
thermal effect (radiative cooling) induces evacuation of the 
flux concentration. The evacuation leads to a downward 
shift of the optical depth scale within the flux concentra- 
tions. The shift is smaller for the smaller features due to 
the late ral radiative heating, w hich inhibits further evac- 
uation ([Venkatakrishnanl Il986[ ) . As a result, the magnetic 
field at r = 1, in the smaller features, is weaker than in the 
case of the feature with more flux. This is in accordance 



with recent numer ical (ICheung et al.l l2007f) and observa- 
tional (jRiiedi et al.lll992l : ISolanki et al.lll996h results. 

During the evacuation, the dowflow velocities reach 
maximum values of 5-10 km/s at r = 1. In the case of 
the largest feature, the downflow extends from the upper 
boundary of the simulation domain and becomes supersonic 
in the lower photosphere. The magnetic features formed 
have diameters of 0.1-0.2 In the case of the biggest fea- 
ture, a supersonic upflow develops approximately 200 s af- 
ter the formation of the flux concentration. The upflow does 
not lead to a complete dispersal of the field, but the feature 
persists until 9 minutes later when it undergoes a fragmen- 
tation. The disappearance of the features in all three cases 
occurs when they meet opposite polarity features, between 
3 and 20 minutes after their formation. 

We also show what happens with the observables when 
the effects of smearing to observational spatial resolution is 
taken into account. An important result is that, in the case 
of small features, Hinode/SP would miss the bright point 
formation and, in some cases, also the high velocity down- 
flows that develop in the process. On the other hand, the 
signatures of the evolution of large features are detectable 
even after the spatial smearing. We show that this case 
can be quantitatively compared with Hinode /SP observa- 
tion (jNagata et al.ll2008l : IShimizu et al.ll2008[ ) and exhibits 
a very similar evolution. This suggests that the magnetic 
field intensification process in the MURaM simulations is a 
faithful description of the process taking place on the Sun. 
Furthermore, our study indicates that the analysis and in- 
terpretation of the observations in terms of the convective 
intensification process is well-founded. 
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